An Extension of the Qz Algorithm beyond the Hessenberg-upper Triangular Pencil
نویسنده
چکیده
Recently, an extension of the class of matrices admitting a Francis type of multishift QR algorithm was proposed by the authors. These so-called condensed matrices admit a storage cost identical to that of the Hessenberg matrix and share all of the properties essential for the development of an effective implicit QR type method. This article continues along this trajectory by discussing the generalized eigenvalue problem. The novelty does not lie in the almost trivial extension of replacing the Hessenberg matrix in the pencil by a condensed matrix, but in the fact that both pencil matrices can be partially of condensed form. Again, the storage cost and crucial features of the Hessenberg-upper triangular pencil are retained, giving rise to an equally viable QZ-like method. The associated implicit algorithm also relies on bulge chasing and exhibits a sort of bulge hopping from one to the other matrix. This article presents the reduction to a condensed pencil form and an extension of the QZ algorithm. Relationships between these new ideas and some known algorithms are also discussed.
منابع مشابه
Performance of the QZ Algorithm in the Presence of Infinite Eigenvalues
The implicitly shifted (bulge chasing) QZ algorithm is the most popular method for solving the generalized eigenvalue problem Av = Bv. This paper explains why the QZ algorithm functions well even in the presence of innnite eigenvalues. The key to rapid convergence of QZ (and QR) algorithms is the eeective transmission of shifts during the bulge chase. In this paper the mechanism of transmission...
متن کاملA Householder-based algorithm for Hessenberg-triangular reduction∗
The QZ algorithm for computing eigenvalues and eigenvectors of a matrix pencil A − λB requires that the matrices first be reduced to Hessenberg-triangular (HT) form. The current method of choice for HT reduction relies entirely on Givens rotations partially accumulated into small dense matrices which are subsequently applied using matrix multiplication routines. A non-vanishing fraction of the ...
متن کاملAlgorithms for Hessenberg-Triangular Reduction of Fiedler Linearization of Matrix Polynomials
Smallto medium-sized polynomial eigenvalue problems can be solved by linearizing the matrix polynomial and solving the resulting generalized eigenvalue problem using the QZ algorithm. The QZ algorithm, in turn, requires an initial reduction of a matrix pair to Hessenberg– triangular form. In this paper, we discuss the design and evaluation of high-performance parallel algorithms and software fo...
متن کاملNumerically Reliable Computation of Characteristic Polynomials
In this paper, we present an algorithm for computing the characteristic polynomial of the pencil (A ? sE). It is shown that after a preliminary reduction of the matrices A and E to, respectively, an upper Hessenberg and an upper triangular matrix, the problem of computing the characteristic polynomial is transformed to the solution of certain triangular systems of linear algebraic equations. We...
متن کاملParallel Reduction of a Block Hessenberg-Triangular Matrix Pair to Hessenberg-Triangular Form—Algorithm Design and Performance Results
The design, implementation and performance of a parallel algorithm for reduction of a matrix pair in block upper Hessenberg-Triangular form (Hr, T ) to upper Hessenberg-triangular form (H, T ) is presented. This reduction is the second stage in a two-stage reduction of a regular matrix pair (A, B) to upper Hessenberg-Triangular from. The desired upper Hessenberg-triangular form is computed usin...
متن کامل